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Abstract 

Recursion relations are presented that allow exact calculation of canonical 
and microcanonical partition functions of degenerate Fermi systems, assuming 
no explicit two-body interactions. Calculations of the level density, sorted by 
angular momentum, are presented for 56 Ni. The issue of treating unbound 
states is also addressed. 
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Statistical models have played a central role in the study of the nucleus since Bohr's 
model of the compound nucleus in 1936 Level densities and occupation probabilities 

are important in determining reaction rates at finite temperature in astrophysical studies 
@ and heavy-ion reactions Q. Level densities also dominate the behavior of low-energy 
neutron scattering ||. A nearly identical set of problems appears in atomic physics when 
modeling highly excited electronic configurations in heavier atoms. There, level densities and 
occupation probabilities play a pivotal role in understanding the opacity of interstellar dust 
which affects the interpretation of a variety of astrophysical measurements |||7j . Mesoscopic 
systems are also beginning to play an increasingly important role in other fields, such as the 
study of quantum dots or atomic clusters. In this letter, we present methods to perform 
exact calculations of both canonical and microcanonical quantities. These methods involve 
the counting of all possible arrangements of fermions amongst independent single-particle 
levels. We apply these methods to the example of counting levels in 56 Ni. Even though 
explicit residual interactions are ignored, the methods allow insight into a variety of problems 
involving finite Fermi systems. We also discuss possible extensions of these pictures that 
would treat residual interactions. 

The calculations presented here are based on recursive algorithms. The algorithms can 
be considered as extensions to techniques used for multi-particle symmetrization which has 
been applied in the context of heavy ion collisions for both Bose and Fermi systems |||| 



and extended to weakly interacting Bose gases [10|. The simplest application considers the 
canonical partition function, for a system of A identical fermions populating a finite set 
of energy levels, e i: at an inverse temperature j3. 
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Here, the A particles are treated as distinguishable. They occupy the states, ji, ■ ■ ■ ]a which 
are eigenstates of H, and V(j) refers to the A\ permutations of j\ ■ ■ ■ ja- 

The permutations may be categorized diagrammatically. As an example, a 9-particle 
permutation is illustrated in Fig. [I], where the n th position in the permutation is found by 



following the segment of the loop that leaves from from the n th site on the diagram. The 
sum over diagrams represents the sum over permutations, and a given closed loop represents 
a permutation cycle. 
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All diagrams can be expressed as a product of cycles. The diagram written in Fig. 1 may 
be expressed as the product, C\C\C±. The sign (— l) n_1 arises from the number of pairwise 
permutations required to create the loop, and is ignored in the case of Bosons. 

One may now derive the recursion relation for the partition function by considering 
summing over diagrams according to the cycle of the loop connecting the upper left point. 
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The factor (A — 1) • • • (A — n+ 1) in the second line of Eq. ([|) may be understood by counting 
the number of ways in which a loop of size n may be connected to one point. 

Eq. (|3|) allows practically instantaneous calculation of the partition function given the 
single-particle energy levels 6j. This approach also leads to a simple expression for two-point 
functions, 
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The four-point function might also be of interest in calculations of residual interactions. 
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It is straight-forward to extend these calculations to arbitrary n-point functions. 

The method is easily extended to include conserved quantities. For instance, conservation 
of the z-component of angular momentum J z can be accommodated by summing over only 
those diagrams which yield a fixed sum, M = Si^i^i- 
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In a similar fashion one may incorporate any other conservation laws, such as parity by 
adding an additional subscript to Z. 

By setting (3 to zero, one may associate Z with a counting of all states regardless of 
energy. By setting (3 = 0, and assigning a discretized energy to each state, one may count 
the net number of states with total energy E. 
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The above relation requires that energy be measured in discrete units. These units can be 
arbitrarily small to make the calculation approach the continuum limit. The quantity N 
represents the micro canonical partition function, and when divided by the energy bin size, 
yields the density of states. 

One may also calculate the microcanonical equivalent of the two and four-point functions. 
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Calculations proceed more quickly when particles and holes are considered independently, 
then convoluted to find the number of states at a fixed E* and M. Considering only states 
with energies above the Fermi level, one first finds N% E * M , the number of ways to arrange a 
particles above the Fermi surface such that they sum to energy E* and angular momentum 
projection M. Then, considering only states below the Fermi surface, one performs the 
analogous calculation for a holes, N% E , M . The number of ways to arrange particles and 
holes subject to conservation of particle number, energy and angular momentum is 

Ne*,M = ^a,E'-E f ,M-M'^a,E*-(E'-E f ),M'- (10) 
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Instead of counting states with a specific projection M = (J z ), one may count Nj, the 
number of multiplets of a given J, by converting the recursion relations for Nm into a relation 
for Nj. 
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By manipulating Eq. (|7]) and Eq.s (|ll|), one may derive a recursion relation for Nj. 
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where the sum over i represents a sum over multiplets of a given j rather than a sum over 
individual states. 

The simplest example to consider is that of equally-spaced non-degenerate levels which 



was addressed by Ericson ||11|| . This example can be related to the problem of finding the 
number of ways an integer can be written as a sum of smaller integers without using the 
same integer twice. The resulting density of states which is displayed in Fig. @ is well 



matched by the Bethe formula H4[L3f . 



exp (27rJgE*/6) 
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where g is the single-particle level density and E* is the excitation energy above the ground 
state. Fig. |3| displays the two-point function, fi = (ajaj) at a fixed excitation energy of 
20 g^ 1 . The resulting phase space filling factor is fairly well described by a Fermi function, 



represented by a solid line in Fig. ^, with the temperature set at T = y^j- 

We now consider a more realistic example where single-particle levels for 56 Ni were gen- 
erated from a Hartree-Fock calculation [[14]]. An infinite external potential was added to 
confine the particles to within a 10 fm sphere, transforming the continuum into a set of 
densely packed discrete states. All levels within 50 MeV of the Fermi surface were included. 
The lowest 20 single-particle energy levels are listed in Table Q. 



TABLES 

TABLE I. Single-particle energy levels of 56 Ni from Hartree-Fock calculations in MeV. 



level 


Proton 


Neutron 


level 


Proton 


Neutron 


°1/ 2 


-32.59 


-42.61 


yy/2 


3.93 


-5.33 


.F3/2 


-24.23 


-34.10 


^5/2 


8.17 


-0.53 


Pl/2 


-22.56 


-32.48 


<^3/2 


9.36 


0.94 


d-5/2 


-15.33 


-25.04 


Sl/2 


9.41 


1.32 


s l/2 


-11.46 


-21.28 


97/2 


12.11 


3.30 


dz/2 


-11.10 


-20.87 


hn/2 


14.37 


5.52 


A/2 


-6.01 


-15.53 


P3/2 


21.42 


7.68 


P3/2 


-1.90 


-11.41 


Pl/2 


14.53 


8.16 


Pl/2 


-0.52 


-9.92 


A/2 


14.95 


8.66 


A/2 


0.80 


-8.67 


A/2 


15.92 


10.00 



The Fermi energy was set at -15.5 MeV for neutrons and -6.0 MeV for protons. The 
energy was discretized in units of 0.25 MeV. One may then couple the resulting level density 
for neutrons with the proton level density to obtain the level density of the composite system. 

Figure [| displays the resulting density of states for three cases. First, the density of 
states is displayed for the case where all single-particle levels are included, even those that 
are unbound. Secondly, results are shown for a calculation where only the fp and g$/2 shells 
are considered, as was implemented in a Monte Carlo shell-model calculation of Nakada and 
Alhassid [15]. Ormand JT^], Langanke [17] and Dean [18] have performed similar studies. 



Continuum levels should not be included in the treatment [HJ. Therefore, the recursive 



sums were modified to subtract the contribution from the continuum in the following manner, 

N A,E = J E^ 1 )^ 1 (E N A - n , E - ei - E N A - n , E ^}j , (14) 

where the primed sums are over eigenstates of the infinite confining spherical well without 
interactions. In the limit that the radius of the confining well approaches infinity, the 
correction should be exact. Figure |] displays the corrected level density for the 56 Ni example. 
The correction only becomes substantial above excitations of 25 MeV. At lower energies the 
aforementioned subset of levels seems sufficient to reproduce the density of states. 

Figure [| displays the number of multiplets within specified energy windows for 56 Ni as 
a function of J. Describing the coupling of angular momentum as a diffusive random walk 



11], one expects the J dependence to behave as 



N E ,j<x (J+~)exp-/?J(J + l), (15) 

where f3 is related to the number of particles and the average j of the single-particle levels 
used to form the angular momentum [|TTJ. As seen in Fig. |5] where f3 was chosen to best fit 
the results, the functional form is accurate except at high J. 

The recursion relations presented here appear quite powerful in calculating statistical 
quantities of excited nuclei with minimal calculational time. The method may easily accom- 
modate an arbitrarily large number of single-particle levels. The 56 Ni calculations required 



approximately 5 minutes of CPU time on a low-end workstation, and most of that time 
was spent performing the coupling of particles to holes described in Eq. QTTD. Enumeration 
of states and their angular momentum is required for nuclear shell-model calculations. In 
current shell-model codes, this is accomplished by combinatorially considering all ways in 
which n particles can be placed into m states. These methods are sufficient as long as a small 
number of single-particle levels are involved. For large nuclei or as one enters or approaches 
the continuum, the number of single-particle levels becomes large and the benefits of the 
recursive method are more apparent. Since the recursive method yields partition functions, 
it is straight-forward to extend it to the calculation of any statistical quantity. 

Although the relations are flexible in handling any model given the single-particle levels, 
two-body (residual) interactions are explicitly ignored. However, they can be included per- 
turbatively. For example, four-point terms can be calculated using Eq.s (||) and (Q). The 
four-point terms can then be used to calculate the contribution of the residual interaction in 
lowest order. It is our hope that the recursive techniques presented here might be expanded 
to include a framework for calculating successively higher order perturbations. 
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FIGURES 




FIG. 1. A sample permutation is represented by the diagram above. The permutation of j\ ■ ■ ■ jg 
in this case is j2,ji, jg, 37, Ja, j&ijs-, jb, jz a s the segment leaving the dot labeled n, is connected to 
j' n in the permutation. This diagram can be expressed in terms of cycles by CiCf C4. 




FIG. 2. The level density is displayed for the case of uniformly spaced levels. It is well ap- 
proximated by the Bethe formula when the excitation energy is much greater than the inverse 
single-particle level spacing, g^ 1 . 
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FIG. 3. The occupation probability is shown for the case of uniformly-spaced levels where the 
energy is fixed at 10 times the single-particle level spacing. 
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FIG. 4. Calculations for the density of states for 56 Ni are displayed for three sets of states. 
Results which consider only fp and the gg/ 2 shell are represented by triangles, including all states 
is represented by squares, while subtracting continuum corrections yields the results represented 
by circles. 




FIG. 5. The number of multiplets of a given J for 56 Ni are shown for the energy windows labeled 
above. The fits to the form shown in Eq. ( |l~5| ) are shown with lines. The form well describes the 
result, except for high J. 



